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Fig.  2.  (a)  Delivery  of  the  MTS  biaxial  tension/torsion  load  frame;  (b)  ASU  graduate  students 

being  trained  by  an  MTS  engineer. 


3.  Approach  and  Accomplishments 


3.1  Optimization  of  specimen  geometry 

Specimen  chosen  fw  testing  in  this  frame  is  a  modified  cruciform  sample  This 
erne, form  spec, men  was  optimized  for  in-plane  fensfon-fension  fesfs  ,o  achieve  a  „  form  s, re 
field  ,n  the  gage  area.  Tubular  specimens  were  fabricated  for  fension/torsion  tests  but 
emp  asis  ,s  being  placed  on  the  cruciform  samples.  The  optimization  process  is  described  next 

3.1.1  Metallic  cruciform  specimens 

The  two  basic  geometries  considered  initially  for  the  cruciform  samples  are  shown  in  Fig.  4. 


(a)  specimen  without  slots 


(b)  specimen  with  slots 


Fig.  4.  FEA  model,  load  and  boundary  conditions 


These  geometries,  which  were  down-selected  from  several  found  in  the  literature  were  studied 

a  rr:  nrr  “e,eiemeni  anaiys,s  <fea>'  wi,h  «*—  -  •» 

and  the  thickness  of  the  gage  section  in  the  cruciform  sample.  Different  corner  fillets  were 
simulated  fo  avoid  local  stress  concentrations.  Various  gage  section  geomeZ  an^oU  n 
each  arm  were  modeled  to  maintain  a  stress  field  as  uniform  as  possible  in  the  gage  ^ 
Due  to  the  symmetry,  only  a  quarter  of  the  specimen  is  analyzed  and  symmetTb^nZ 
conditions  are  applied.  The  ratio  of  applied  loads  along  the  y  (vertical)  and  x  (horizontal)  axes  Z 
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Py/Px  -  1.0.  Figure  5  shows  contours  of  resultant  displacement  in  the 
simple  and  slotted  cruciform  specimens  obtained  from  FEA. 


gage  area  of  both  the 


(a)  specimen  without  slots  (b)  specimen  with  slots 

Fig.  5.  Displacement  distribution  from  FEA. 


igure  shows  the  Von  Mises  stress  distribution  for  both  geometries.  It  can  be  seen  that  the 

and1  the  ni'T^f  °l!  Sh°WS  concentrations  at  the  corner  between  the  gage  area 

and  the  fillet  wh.le  the  stress  distribution  is  uniform  for  the  specimen  with  the  slotted  arms 


Fig.  6.  Von  Mises  stress  distribution  from  FEA. 


Figures  7  and  8  show  the  maximum  principal  stress  and  in-plane  shear  stress  distributions 

rry-  !he  resu,,s  indica,e  ,ha* ,he  «*""  «,»,  ; 

generate  an  almost  perfect  two  dimensional  stress  state  within  the  gage  area;  therefore  that 
geometry  was  chosen  for  metallic  samples. 


(a)  specimen  without  slots  (b)  specimen  with  slots 

Fig.  7.  Maximum  principal  stress  distribution  from  FEA. 
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(a)  specimen  without  slots  (b)  specimen  with  slots 

Fig.  8.  In-plane  shear  stress  distribution  from  FEA. 


The  final  sample  geometry  and  an  example  of  a  machined 


specimen  are  shown  in  Fig.  9. 
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Fig.  9.  (a)  Final  dimensions  (in  inches)  of  the  cruciform  sample;  (b)  Machined 


specimen. 


A  modified  cruciform  specimen  with  a  central  hole  (Fig.  10)  was  also  modeled  under  identical 
toads  and  used  to  study  fatigue  crack  initiation  at  the  inner  circumference.  The  FEA  shows  a 
Mises  stress  concentration  around  the  inner  circumference  that  decreases  away  from  the  hoTe 


Fig.  10.  von  Mises  stress  distribution  of 


a  modified  cruciform  specimen  with  a  central  hole. 
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The  basic  geometry  shown  above  was  further  modified  to  create  a  specimen  that  could  be  used 
study  correlations  between  microstructure  and  fatigue  crack  nucleation  under  biaxial  loadino 

“  ,  3 '  aS  We"  as  reSUl,S  ,rom  FEA  °f  lhe  specimen  to  verify  thZ  he 

concentration  around  the  centra,  ho,e  is  indeed  uniform  despite  the  reduced  sample  lize 


3.1.2  Cruciform  specimens  for  composite  materials 

loaZotn  StmC,UreS  are  0f,en  SUb,eC,ed  '°  comP,icated  service  loading  conditions  These 

load,™  ™  are, Very  rare'y  UniaXia';  ,herefore'  ,w°-  °r  three-dimensional  static  and  cyclic 
loading  ,s  necessary  for  realistic  testing.  In  order  to  insure  reliable  operation  and  efficient  destan 

loaT^mP,°nen>S'  "  'S  Cri,iCal  10  charac,erize  the  material  behavior  under  complex  multiaxL 

i  i  z  be;sed  ,or  s,ruc,urai  hea,m  “*• 

composite  sZdmens.  ^  3  56,168  °'  '6S,S  Wera 

™i,T"Z',Ca,Cl,0rm  sp80lmen  also  needad  K>  be  designed  in  this  case.  The  design  criteria 

shearltra  nsTnZ'oao0"  '"rf  ^  r69i°n  Uniform  biaXial  S,rainS'  (li)  mlnlml2ation  of  the 
section  '  <l")  min,miza,ion  °<  ,ha  stress  concentrations  outside  of  the  gage 

ZrZ« ‘i9a,e  k6  °'  9e0rne,rV  °n  ,he  s,ress  concentration  of  the  specimen  geometric 

thicZsTnfm  35  '  'he  rat"US  °f  *he  “mer  f"el  at  ,he  ibtersection  of  the  arms  (ii)  the 
ckness  of  the  gage  area  and  (ill)  the  geometry  of  the  gage  area  were  considered  Some  of  the 

e  fferent  geometries  that  were  fabricated  and  tested  are  shown  in  Figures  12  3  and  H  A I  of 
them  included  a  central  hole  to  induce  crack  nucleation  at  a  controlled  location. 
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y  (b) 

Fig.  12.  Uniform  thickness  sample,  (a)  Geometry  (in  inches);  (b)  Fabricated  specimen. 
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Fig.  13.  Variable  thickness  sample. 


(b) 

(a)  Geometry  (in  inches);  (b)  Fabricated  specimen. 
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Fig.  14.  Variable  thickness  sample,  (a)  Geometry  (in  inches);  (b)  Fabricated 


specimen. 


The  gage  area  of  the  specimen  shown  in  Fig.  12  is  the  same  thickness  as  the  arms  and  it  has  a 
atea  Z  ‘  ’  °'  inte^  of  the  arms,  while  the  thickness  of  the  gage 

,025J“  an<)  M  S  1/3  °fth8  lhlcknesso,,he  arms.  and  15.24  mm  (0.6")  and  6.35  mm 

(0.25  )  radii  for  the  corner  fillets,  respectively.  Cracks  should  initiate  at  the  central  hole. 
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3.2  Tubular  specimens 

Tubular  specimens  were  designed  for  tension/torsion  as  shown  in  Fig.  15.  This  loading  has  the 
advantage  of  providing  variable  principal  directions  for  non-proportional  loading.  Figure  16 

shows  a  specimen  subjected  to  a  tension/torsion  cyclic  loading.  It  can  be  seen  that  failure 
happened  at  the  gage  section,  as  expected. 
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Fig.  15.  Tubular  specimens,  (a)  Geometry  (in  inches);  (b)  Fabricated  samples. 


Fig.  16.  Broken  tubular  specimen  under  tension/torsion  cyclic  loading 

4.  Applications  of  the  Biaxial  Frame 

In  this  section,  some  applications  of  the  biaxial  frame  will  be  reported  for  material  and  structural 
damage  detection  and  residual  useful  life  estimation. 


4.1  Fatigue  crack  monitoring  in  metallic  samples 

The  modified  metallic  cruciform  specimen  with  central  hole  was  tested  under  cyclic  loading  with 
a  load  ratio  R  =  0.1  and  frequency  of  10  Hz.  The  maximum  applied  load  of  32  kN  along  both  the 
x  and  y  axes  was  based  on  the  results  of  the  nonlinear  FE  simulation,  so  that  the  von  Mises 
stress  around  the  hole  was  slightly  lower  than  the  yield  strength  of  Al  2024,  which  was 
measured  from  a  uniaxial  test  to  be  325  MPa.  The  experimental  test  setup  for  the  cruciform 
specimen,  including  the  biaxial  frame,  the  acoustic  emission  (AE)  system,  the  data  acquisition 
system  (DAQ)  and  the  CCD  camera,  is  shown  in  Fig.  17. 
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Fig.  17.  Experimental  setup  for  damage  monitoring  in  a  metallic  cruciform  specim. 
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Figure  1 8  shows  the  placement  of  the  AE  and  piezoelectric  (PZT)  transducers.  In  this  setup  the 
AE  sensors  were  used  for  crack  monitoring  and  the  PZT  sensors  for  damage  diagnosis  and 
prognosis.  Figure  19  shows  the  CCD  camera  used  for  monitoring  crack  propagation. 


Fig.  18.  Placement  of  AE  and  PZT  sensors  in  a  metallic  cruciform  sample. 
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Fig.  19.  The  installation  of  CCD  camera 


Figure  20a  shows  a  1  mm  long  crack  after  1598951  cycles.  The  AE  amplitude  then  suddenly 
increased  to  66  dB,  as  shown  in  Fig.  20  (b),  which  was  probably  due  to  rapid  crack  growth 
through  the  thickness  of  the  specimen,  as  the  initial  crack  probably  had  a  thumbnail  geometry 


At  1610561  cycles,  the  crack  reached  a  4  mm  length  and  the  AE  amplitude  increased  to  76  dB 
as  shown  in  Fig.  21  (a)  and  (b).  It  should  be  noted  that  the  crack  initiation  could  occur  at  any 
place  around  the  inner  circumference  of  the  hole  and  the  crack  could  grow  along  either  the  x  or 
the  y  direction  based  on  the  nonlinear  FEA  simulation  shown  in  Fig.  11.  The  crack  growth 
direction  showed  in  Figures  20  and  21  is  close  the  y  direction,  probably  because  it  was  parallel 
to  the  rolling  direction  of  the  2024  Al  plate  used  to  fabricate  the  sample,  and  the  rolling  direction 
tends  to  be  weaker  with  respect  to  crack  propagation. 
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Fig.  21.  (a)  Fatigue  crack  4  mm  long  observed  at  1610561  cycles;  (b)  AE  amplitude  of  76  dB. 


4.2  Structural  health  monitoring  (SHM)  on  advanced  composites 

Although  composite  materials  are  being  used  in  various  structural  systems  such  as  aircrafts  and 
au  omobiles,  SHM  and  damage  prognosis  of  composite  materials  are  still  emerging  research 
areas.  This  section  will  report  on  the  use  of  the  biaxial  frame  for  SHM  of  woven  composites 


4.2.1  Composite  materials  used  in  the  test 

Multiple  carbon  fiber  polymer  composite  test  coupons  were  fabricated  in  the  AIMS  center  at 

Un,VerS,ty-  Woven  fiber  was  Purchased  from  Fiber  Glast  Inc.  and  the  matrix  used 
was  HEXION  EPON  863  and  EPI-CURE  3290.  A  hot  press  was  used  to  cure  composite  plates 
and  the  shapes  of  different  coupons  were  cut  using  CNC  machines.  Different  stacking 
sequences  were  also  investigated  to  see  how  it  affected  the  fatigue  life  of  the  coupon. 


4.2.2  Uniaxial  testing 

Uniaxial  tensile  tests  were  completed  to  determine  the  tensile  strength  of  the  materials  used  in 
a  E'9ht'ply’  Puasi-isotropic  woven  composite  sample  was  used  in  the  test  with 

[0  /45  /-45  /90  ]s  stacking  sequence.  The  monotonic  tensile  test  followed  the  ASTM  standard 
D  638.  The  geometry  of  the  test  coupon  and  examples  of  pristine  and  damaged  samples  are 
shown  in  Fig.  22.  The  measured  strength  of  the  composite  coupons  was  used  to  estimate  the 
cyclic  load  for  the  biaxial  tension-tension  fatigue  test. 
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Fig.  22.  (a)  Geometry  (in  inches)  of  the  composite  dog  bone  sample,  (b)  Uniaxial  loading  test 

using  composite  dog  bone  samples. 


4.2.3  Biaxial  Test 

The  biaxial  fatigue  tests  for  composite  cruciform  specimens  were  conducted  to  study  the 
pe  ormance  of  composite  structures  under  complex  loading  conditions.  In  this  work,  multiple 
sensing  techniques,  such  as  AE,  PZT  active  and  passive  sensing  and  strain  gages  were  used  to 

omtor  the  degradation  process  of  the  composite  material.  The  experimental  setup  for  the 
biaxial  test  is  shown  in  Fig.  23. 


Fig.  23.  Testing  setup  for  composite  cruciform  samples  under  biaxial  fatigue  load. 


™  STSOr  S'9nalS  we,e  collected  a9a'nst  three  different  types  of  actuation  inputs: 
narrow  band  burst  tnpuf,  wide  band  chirp  input  and  pseudo  random  pulse  trains.  The  signal 
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eatures  from  various  active  and  passive  methods  were  being  evaluated  in  the  context  of 
different  metrics  to  eventually  develop  a  high  fidelity  and  robust  prognostic  algorithm.  A  time 
series  prognostics  algorithm  based  on  a  multivariate  Gaussian  process  was  used  to  map  the 
input  sensor  signal  features  to  a  time  varying  damage  index.  In  addition,  the  algorithm  was  used 
o  orecast  the  future  damage  index  and  the  corresponding  residual  useful  life  estimate  (RULE). 

e  prognostics  framework  will  be  verified  by  experiments  with  different  biaxial  loadings.  The 
experimental  excitation  and  sensing  signals  are  shown  in  Figures  24  to  26. 


Fig.  24.  Chirped  excitation  signals  and  related  sensing  signals. 


Fig.  25.  Burst  excitation  signals  and  related  sensing  signals. 
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Excitation  Pseudo  Rar  dom  Pulse  Trans  Sensing  Signals  of  Sensor  t 


Fig.  26.  Pseudo  random  pulse  trains  and  related  sensing  signals 


Acoustic  emission  techniques  were  also  used  during  this  test.  Four  sensors  were  mounted  on 
the  coupon  and  the  counts  associated  to  acoustic  events  indicated  a  clear  trend  in  the 
degradation  process.  The  AE  counts  are  a  good  choice  at  later  stages  of  fatigue  life  because  a 
large  amount  of  high-energy  acoustic  bursts  were  released  as  damage  occurred.  In  Fig.  27,  the 
number  of  counts,  amplitude  and  time  are  shown.  Note  the  steady  increase  of  the  amplitude  up 
to  92  dB  indicating  the  growth  of  damage. 


Fig.  27.  Acoustic  emission  Hits  vs.  Time  (sec)  and  Amplitude  (dB)  during  testing. 


4.3  SHM  on  a  metallic  specimen 

In  this  application,  a  correlation  analysis  approach  based  on  non-parametric  system 
identification  is  used  to  estimate  the  current  damage  state.  This  damage  state  is  estimated  by 
correlating  the  real-time  sensor  signals  obtained  from  two  different  strain  gages  placed  at  two 
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different  places  on  the  structure/specimen.  Once  the  damage  state  (in  the  form  of  a  damage 
index)  is  estimated,  it  is  fed  to  an  off-line  predictive  model  based  on  a  Gaussian  process  to 
forecast  the  future  damage  state  and  the  residual  useful  life.  The  real-time  algorithm  was 
validated  on  an  AI-6061  cruciform  specimen  undergoing  biaxial  fatigue  loading.  The  details  on 
the  application  of  this  procedure  are  provided  in  the  paper  included  in  the  appendix. 

5.  Outreach 

5.1  Local  schools 

The  Kyrene  Middle  School  Prep.  Program  (8th  grade),  AZ,  visited  the  Biaxial  Testing  Facility  on 
4/30/09.  These  students  were  given  a  laboratory  tour  and  some  basic  instruction  (visual)  on 
multi-axial  loading  using  the  biaxial  test  frame.  A  picture  is  shown  in  Fig.  28. 


Fig.  28.  Kyrene  Middle  School,  AZ,  students  visit  the  biaxial  testing  facility. 

5.2  AFOSR  sponsored  research 

The  program  managers  of  the  SHM  MURI  project,  Dr.  Victor  Giurgiutiu  and  Dr.  David  Stargel 
from  AFOSR,  and  members  of  the  MURI  Advisory  Board  visited  the  Biaxial  Testing  facility  on 
5/19/09.  They  were  given  a  demonstration  of  the  use  of  the  frame  on  SHM  of  metallic  materials 
and  were  shown  the  corresponding  testing  setups,  as  shown  in  Fig.  29. 


Fig.  29.  MURI  program  managers  and  advisory  board  members  visit  the  biaxial  testing  facility. 
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Fatigue  Testing."  Proceedings  of  the  7th  International  Conference  on  Fatigue  Damage  of 
Structural  Materials,  Sept.  2008. 

7.  Summary 

■  Optimized  cruciform  specimens  for  metallic  materials 

«  Optimized  cruciform  specimens  for  advanced  composite  materials 

■  Designed  and  fabricated  miniature  cruciform  samples  to  study  crack  nucleation 

■  Designed  tubular  metallic  specimen  for  tension/torsion  testing 

■  Conducted  crack  initiation  test  via  Al  2024  cruciform  specimens 

■  Conducted  SHM  validation  test  via  Al  2024  cruciform  specimens 

■  Conducted  SHM  validation  test  via  woven  composites  cruciform  specimens 

8.  Appendix 
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On-line  Structural  Health  Monitoring  and  Prognosis  of 
a  Biaxial  Cruciform  Specimen 


Subhasish  Mohanty  *  Aditi  Chattopadhyay  t  Jun  Wei  ^ 

Pedro  Peralta5 

Mechanical  and  Aerospace  Engineering,  Arizona  State  University,  Tempo,  AZ ,  85287,  USA 


The  current  research  of  on-line  damage  state  estimation  techniques  offers  adaptive  dam¬ 
age  state  prediction  and  residual  useful  life  assessment.  The  real-time  damage  state  infor¬ 
mation  from  an  on-line  state  estimation  model  can  be  regularly  fed  to  a  predictive  model  to 
update  the  residual  useful  life  estimat  ion  in  the  event  of  a  changing  situation.  The  present 
paper  discusses  the  use  of  an  integrated  prognosis  model,  which  combines  an  on-line  state 
estimation  model  with  an  off-line  predictive  model  to  adaptively  estimate  the  residual  use¬ 
ful  life  of  an  AI-6061  cruciform  specimen  under  biaxial  loading.  The  overall  fatigue  loading 
history  is  assumed  to  be  a  slow  time  scale  process  compared  to  the  time  scale  at  which, 
the  sensor  signals  are  acquired  for  on-line  state  estimation.  The  fast  scale  on-line  model  is 
based  on  a  non-para  me  trie  system  identification  approach  such  as  correlation  analysis.  A 
new  damage  index  equivalent  to  quantitative  damage  state  information  at  any  particular 
fatigue  cycle,  is  proposed.  The  on-line  model  regularly  estimates  the  current  damage  state 
of  the  structure  based  on  passive  strain  gauge  signals.  These  damage  states  information  is 
regularly  fed  to  the  slow  scale  off-line  predictive  model  as  it  becomes  available.  The  off-line 
predictive  model  is  a  probabilistic  nonlinear  regression  model,  which  is  based  on  Bayesian 
statistics  based  Gaussian  process  approach.  The  off-line  module  adaptively  updates  the 
model  parameters  and  recursively  predicts  the  future  states  to  provide  residual  useful  life 
estimate. 


N  omenclat  ure 


Fast  scale  variables 
m  Fast  scale  lag  numbers 

At  (—  1///)  Fast  scale  time  interval 

T  Fast  scale  observation  time 

M  No.  of  fast  scale  observation  samples 

u(m)  Fast  scale  input  observation  at  lag  m 

v{m)  Fast  scale  output  observation  at  lag  m 

v{rn)  Fast  scale  noise  observation  at  lag  m 

luv{nh)  Fast  scale  cross  correlation  coefficient  at  lag  m 
Slow  scale  variables 

n  Slow  scale  damage  level  number  (or  damage  instances) 

(Note:  In  general  n  is  not  same  as  number  of  fatigue  cycles) 

*  Number  of  damage  instances  after  last  fast  scale  data  available 

n  After  last  fast  scale  data  available, the  number  of  damage  instances  to,  damage  state  become  critical 

N  Total  no.  of  fatigue  cycles 

N0  No.  of  fatigue  cycles  elapsed,  before  the  first  on-line  sensor  data  available 
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AN  No.  of  fatigue  cycles  per  increment  of  each  damage  instant 

Slow  scale  material  scatter  related  noise  hyperparameter  at  nth  fatigue  cycles 
#£  Slow  scale  process  hyperparameter  at  nth  damage  instances 

0™  Slow  scale  input  weighing  hyperparameter  at  nth  damage  instances 

0bn  Slow  scale  input  bias  hyperparameter  at  nth  damage  instances 

Slow  scale  input  vector  at  nth  damage  instances  with  n  =  1, 2 . . . ,  n, . . .  y  Ar 
Slow  scale  output  damage  index  at  nth  fatigue  cycles  with  n  —  1. 2 . . . ,  n, . . , .  N 
Kn  Slow  scale  kernel  matrix  at  nih  fatigue  cycles  with  n  =  1, 2 . . . ,  n» . . , .  N 

/c(xj,  Xj)  Slow  scale  kernel  function  at  nth  fatigue  cycles  with  i,  j  ==  1. 2 . . .  . ,  n, .. . . .  Ar 


L  Introduction 

Real-time  health  monitoring  and  prognostics  is  emerging  at  the  forefront  of  Condition  based  Maintenance 
(CBM)  of  critical  structural  systems  giving  rise  to  the  term  Prognostic  Health  Management  (PHM).  Whether 
it  is  a  newly  acquired  or  an  aging  fleet  aircraft,  the  structural  life  ceiling  of  the  fleet  aircraft  are  defined  from 
three  distinct  approaches:  safe-life,  fail-safe,  and  dam  age- tolerant  approaches.  A  detailed  review  of  these 
approach  and  their  applications  to  US  Navy’s  P-3C  aircraft  is  presented  by  Iyer,  ct.  ah1  In  the  case  of  safe-life 
approach  service  or  retirement  life  of  a  component  is  defined  by  crack  initiation  time  derived  from  a  full-scale 
component,  or  element  fatigue  test.  In  practice,  however,  the  component  is  retired  before  the  formation  of 
a  fatigue  crack  by  using  safety  factors  on  calculated  crack  initiation  time  because  of  inherent  variability 
in  both  static  and  fatigue  material  properties  as  well  as  assumptions  made  in  analytical  models  used  to 
calculate  crack  initiation  time.  The  fail-safe  approach  assumes  an  initial  damage  due  to  manufacturing  and 
its  subsequent  growth  during  service  to  detectable  crack  sizes.  Service  life  in  fail-safe  structures  can  thus  be 
defined  as  the  time  to  a  survivable  detectable  damage.  The  fail-safe  approach  requires  inspections  as  part 
of  the  maintenance  program,  and  these  inspections  are  also  specifically  geared  to  find  damages  in  identified 
critical  components  of  the  aircraft.  In  the  case  of  damage-tolerant  approach,  it  assumes  initial  defects, 
however  small  they  may  be  in  critical  structural  elements,  which  will  eventually  grow  in  service  to  large  crack 
sizes.  Service  life  is  estimated  through  rigorous  crack  growth  analysis  both  deterministic  and  stochastic.  A 
majority  of  the  deterministic  crack  growth  analysis  models  are  based  on  the  approach  followed  in  FASTRAN2 
or  AFGROW3  type  crack  growth  model.  It  is  to  be  noted  that  these  models  are  based  on  some  empirical 
parameters,  which  have  to  be  tuned  for  a  particular  application.  In  case  of  stochastic  life  estimation,  two 
different  approaches,  are  currently  being  used  in  industry,  and  are  these  techniques  are  continuously  being 
improved  by  the  stochastic  life  modeling  research  community.  In  the  first  technique  service  life  is  based 
on  the  basic  principles  of  Miner’s  damage  accumulation  criteria,  and  in  the  second  technique  the  stochastic 
crack  growth  curves  are  obtained  first,  followed  by  estimation  of  the  corresponding  stochastic  service  life. 
Wu  and  Ni4  presented  a  stochastic  crack  growth  model  which  can  be  used  for  stochastic  life  estimation 
of  structures.  In  the  case  of  damage  accumulation  type  models,  the  recent  work  by  Liu  and  Mahadeven5 
on  nonlinear  fatigue  damage  accumulation  rule  and  a  stochastic  S-N  curve  representation  technique  for 
predicting  stochastic  fatigue  life  under  variable  amplitude  loading  is  noteworthy.  In  the  above  mentioned 
deterministic  or  stochastic  approaches  the  damage  tolerance  and  fatigue  life  predictions  are  obtained  based 
on  assumed  structural  flaws  or  based  on  previous  coupon  test  results  regardless  of  whether  the  assumed 
structural  flaws  actually  occur  in  service  (in.,  prognosis  is  made  before  diagnosis).  In  addition  in  a  real  life 
scenario,  changing  loading  conditions  and  other  environmental  conditions  such  as  change  in  humidity  and 
temper atme,  leads  to  different  Crack  growth  law  parameters,  compared  to  those  originally  found  from  coupon 
testing.  Consequently,  a  large  degree  of  conservatism  is  incorporated  into  structural  designs  due  to  these 
uncertainties.  The  current  research  in  the  area  of  on-line6"8  damage  state  estimation,  or  Structural  Health 
Monitoring  (SHM)  techniques  offers  methodologies  for  adaptive  damage  state  prediction  and  residual  useful 
life  assessment.  The  on-line  state  estimation  model  can  be  either  supervised  or  baseline-free  unsupervised 
model.  In  a  supervised  approach,7,8  the  model  has  to  be  trained  for  known  damage  cases  based  on  previously 
conducted  tests.  The  supervised  model  is  difficult  to  implement  in  real-life  condition.  This  is  because  of 
the  changing  boundary  conditions  and  noise  levels  (due  to  electrical  connections),  which  may  not  necessarily 
be  similar  to  the  previously  obtained  data  (those  used  for  training  purpose).  A  detail  review  on  different 
SHM  approaches  are  presented  by  Farrar,  et.  al.9?1°  For  real-time  damage  monitoring  and  prognosis,  there 
is  also  a  need  for  an  effective  predictive  model  to  forecasts  the  future  state  and  the  remaining  life  of  the 
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structure.  The  real-time  damage  state  information  from  the  on-line  state  estimation  model  can  be  regularly 
fed  to  the  predictive  model  to  update  the  residual  useful  life  estimation  in  the  event  of  a  changing  situation 
The  bearing  failure  adaptive  predictive  model  presented  by  Billington,  et.al.12  is  among  the  few  earliest 
reported  work  on  real-time  adaptive  predictive  models.  However,  the  model  is  based  on  a  linear  covariance 
structure,  which  may  not  always  be  suitable  for  learning  nonlinear  damage  growth  dynamics.  The  present 
paper  proposes  a  recursive  Gaussian  process  predictive  model,  in  which  the  model  parameters  arc  adaptively 
updated  to  predict  the  future  states  and  residual  useful  life  estimate.  Unlike  the  covariance  based  structure 
he  Gaussian  process13-15  uses  the  kernel  function.  The  kernel  function  is  a  multi-dimensional  function  and 
transfer  the  linearly  inseparable  information  first  to  a  high  dimensional  feature  space,  where  the  information 
can  be  linearly  separable.  In  the  present  paper  a  non-parametric  system  identification  based  correlation 
analysis  approach  is  used  to  estimate  the  current  damage  state.  The  current  damage  state  is  estimated 
y  correlating  the  real-time  sensor  signals  obtained  from  two  different  strain  gauges  placed  at  two  different 
places  on  the  structure.  Once  the  damage  state  (in  the  form  of  a  damage  index)  is  estimated,  it  is  fed  to 
the  Gaussian  process  off-line  predictive  model  to  forecast  the  future  damage  state  and  residual  useful  life 
estimate  (RULE).  The  real-time  algorithm  is  validated  on  a  Al-6061  cruciform  specimen  undergoing  biaxial 


II.  Integrated  Prognosis  Model 

The  integrated  prognosis  model  is  a  real-time  prognosis  model,  that  work  in  conjunction  with  real-time 
sensor  signal  measurements.  As  shown  in  Fig.  1,  the  integrated  prognosis  architecture  has  two  distinct 
sub  modules,  the  on-line  state  estimator  and  the  off-line  state  predictor.  The  on-line  state  estimator  infers 
the  current  state  of  the  structure  from  real-time  sensor  measurements.  Once  the  current  on-line  state 
information  becomes  available,  it  is  fed  to  the  off-line  predictive  model  to  predict  the  future  states  and 
compute  the  corresponding  residual  useful  life  estimate.  It  is  noted  that  the  overall  fatigue  loading  history  is 
assumed  to  be  a  slow  time  scale  process6  compared  to  the  time  scale  at  which  the  sensor  signals  are  acquired 
at  discrete  intervals.  Based  on  this  assumption,  from  this  point  onwards  the  sensor  signal  acquiring  process 
and  corresponding  current  cycle  damage  state  (or  damage  index)  estimation  process  is  denoted  a  fast  scale 
process,  whereas  the  overall  slower  fatigue  process  is  denoted  a  slow  scale  process.  It  is  noted  that  the 
fast  scale  sensor  measurements  are  performed  at  discrete  slow  scale  intervals.  Also,  it  is  noted  that  at  the 
individual  slow  scale  instances,  where  the  fast  scale  sensor  measurements  are  taken,  the  damage  state  of  the 
structure  is  assumed  to  remain  unchanged. 
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Figure  1.  Schematic  of  integrated  prognosis  model. 
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A.  On-line  fast-scale  damage  state  estimation 

1.  Fast-scale  transfer  function 

At  any  given  nth  slow  scale  fatigue  cycle  the  fast  scale  output  sensor  measurement  can  be  mapped  with  the 
fast  scale  input  sensor  measurement  over  fast  scale  unmeasured  noise  v.  A  typical  nth  fatigue  cycle  block 
diagram  that  maps  the  fast  scale  input  measurement  with  fast  scale  output  measurement  is  shown  in  Fig. 

2.  The  fast  scale  Z-domain  transfer  function  P(z)  between  input  u  and  output  y  at  nth  damage  level  can  be 
represented  as: 


Unmeasured 
electrical  noise 


v(m) 


Measured  input 

Degrading  structure 

signal 

at  n111  damage  level 

7 

u{m) 

- *v 

with  plant:  P„(z ) 

Measured 


y(m) 


Figure  2.  Block  diagram  for  fast  scale  transfer  function.  The  transfer  function  is  an  instantaneous  represen¬ 
tation  of  the  time  degrading  structure  at  any  typical  damage  level.  However  as  the  damage  grow  the  transfer 
function  also  changes  leading  to  a  time  variant  approach  for  system  identification. 


v(t)  =  P(z)u(t)  +  v(t) 

—  {bu  +  b\z  1  4  62-2  2  +  •  •  •  4-  b^jz~M)u(t)  4-  v(t)  (1) 


where  z~m;m  =  0, 1, ..  .  M  arc  the  backspace  operators  of  the  pulse  transfer  function  P(z)  and  bTn\m  = 
0,1,... M  are  the  finite  impulse  response  (FIR)  coefficients.  Then  Eq.  (2)  can  be  rewritten  in  the  discrete 
domain  as 

v(t)  =  b0u(t)  4-  b\u(t  -  1)  4-  b2u(t  -  2)  +  . . .  4-  bMu(t  -  M)  +  v(t)  (2) 


2.  Damage  index 

The  slow  scale  damage  index  a(n),  at  nth  fatigue  cycle,  is  the  representative  damage  state  inferred  from  nth 
fatigue  cycle  fast  scale  sensor  measurements.  The  damage  index  can  be  derived  by  finding  the  mth  lagged 
output  v(t-\-m)  from  Eq.  (2)  and  pre  multiplying  the  input  u(t),  obtaining 

u(t)v{t-\-m)  —  bQu(t)u(t+m)+biu(t)u(t  —  \+rn)+b2u(t)u(t-2+m)  +  . .  .+bMu(t)u(t— M +m)+u(t)i'(t+m) 
.  .  (3) 

Applying  expectation  operator  to  both  sides  of  Eq.  (3)  and  assuming  independence  between  noise  and  the 
input  signal,  the  mth  lagged  cross-correlation  coefficients  can  be  expressed  as 

7 uv(m)  =  b0yuv(m)  +  bnuv(m  -  1)  +  b2yuv(m  -  2)  +  . . .  +  bM^uv(m  -  M)  ;  m  =  0, 1 . . . ,  M  (4) 

With  known  input  (u)  and  output  (v)  time  series,  the  mth  lagged  cross-correlation  coefficients  7uv(^i)  and 
auto-correlation  coefficients  7u(m),  the  FIR  coefficients  6m;m  =  0,1,..  .M  can  be  estimated  and  hence 
the  nth  fatigue  cycle  pulse  transfer  function.  To  estimate  M  -F  1  FIR  coefficients  we  need  to  solve  M  +  1 
algebraic  equations  given  by  Eq.  (4).  However,  solving  M  +  1  algebraic  equations  involves  inverting  a 
(M  +  1)  x  (M  4  1)  autocorrelation  coefficient  matrix,  which  becomes  computationally  expensive  in  the 
context  of  real-time  applications.  To  circumvent  this  problem  rather  than  directly  estimating  the  transfer 
function  and  the  damage  states  (such  as  nth  fatigue  cycle  zeros  and  poles),  the  damage  state  equivalent 
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damage  index  can  be  estimated.  Based  on  the  cross-correlation  coefficients  %v(m)  a  new  damage  index  is 
formulated,  which  is  expressed  as  follows: 


sr-vm—iVf 

2-^ni—0 


1 


E 


•(7gyM-7gyM)2 

ZZo  WyM)2 


;-n - 


1,  2 . . .  ,AT/AN 


(5) 


where  7^v(m)  represents  fatigue  cycle  fast  scale  cross-correlation  coefficients,  and  7^  (m)  represents  0th‘ 
(or  reference)  fatigue  cycle  fast  scale  cross-correlation  coefficients.  It  is  noted  that  if 
max  {ju(m)  or  7 v{m)  —  1:  m  —  0, 1,  . . .  M},  then  either  the  sensor  is  damaged  or  is  debonded,  pro¬ 
ducing  only  noise. 


B.  Off-line  slow-scale  damage  state  prediction  and  residual  useful  life  estimation 

1.  Predicting  in  a  Bayesian  framework 

The  goal  of  a  probabilistic  Bayesian13'15  forecasting  approach  is  to  compute  the  posterior  distribution11  of 
a  future  damage  state  or  damage  index  n7l+i  ,  i.e.,  to  determine  the  probability  distribution  of  the  random 
damage  index  an+\  given  a  random  test  input  and  a  set  of  n  training  data  points  described  a s  D  = 
the  Bayesian  framework  the  predictive  distribution  with  mean  and  variance  can  be  found 

by  conditioning  the  damage  indices  a±)  a2, . that  are  affected  by  the  corresponding  random  inputs 

xi,  x2r..,*.,xu,Xn+i-  A  prior  over  the  space  of  possible  functions  to  model  the  random  damage  index  as 
/(a | a)  ,  can  be  defined  where  a  are  some  hyperparameters  that  can  account  for  random  load  sequence  effect 
in  the  form  of  curve  fitting.  Also,  a  prior  over  the  noise  f(d |/?),  can  be  defined  where  d  is  some  appropriate 
noise  vector  that  arises  due  to  scatter  in  material  micro  structure  and  0  is  another  set  of  hyperparameters 
used  to  mode!  the  uncertainty  due  to  scatter.  Now  if  the  hyperparameters  a  and  0  are  given,  the  conditional 
probability13”15  can  be  expressed  as 

/  (an+il  {Xi=i,...,„,  Q',/3})  =  J (an+i|  {xi=i,...in,  a,  })  f(a\a)f(tf\j3)dad-d 

(6) 

where  a  and  d  denotes  the  underlying  function  which  respectively  corresponds  to  damage  index  and  noise 
due  to  scatter.  Since  ai,  02,......^.  and  an+i  are  conditioned  random  variables  in  the  observed  set  of  damage 

indices,  the  conditional  distribution  of  an+1  can  be  written  as  follows 


2 .  Predicting  with  a  Gaussian  process 


To  evaluate  Eq.  (7)  it  is  necessary  to  evaluate  the  integral  given  in  Eq.  (6).  However,  in  general,  Eq.  (6)  is 
complicated  to  evaluate.  The  standard  approach  to  evaluate  the  integral  in  Eq.  (6)  is  by  a  method  called 
evidence  maximization16  or  by  numerically  integrating  by  Monte  Carlo  simulation.17  However,  assuming  the 
underlying  damage  index  function  follows  a  Gaussian  distribution,  the  exact  analytical  farm 

Of  Eq.  (7)  is  as  follow 

/  (an+l  l  v^n+l)  =  (2'7r)n+172det(K.7J.+1)1^2  2^an’fl  ”  ^  Kn+i  (an+I  “  M)^ 

(8) 


where  a  is  the  function  mean  and  Kn  is  a  n  x  n  kernel  matrix.  The  individual  elements  knmn  of  the  kernel 
matrix  K„  can  be  found  from  a  parameterized  kernel  function  that  will  be  described  in  the  next  section. 
Assuming  zero  mean  function  distribution  Eq.  (8)  can  be  written  as 


(27r)det(K„+1) 


r,)  exp(- 


2(7? 

«n+l 
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where  au+ j  is  the  one-step  ahead  predicted  mean  at  slow  scale  damage  level  n  +  1  and  is  given  by 

an+i  =kTK“Ian;  kx  =  A:(xn+i,xi)feIi2,...„  (10) 

whereas  <?ln+1  is  the  one-step  ahead  predicted  variance  at  slow  scale  damage  level  n  +  1  and  is  given  by 

=  fc-k^IQ'k;  fci  =  fc(xn+1,Xi)i=i(2,..,»;  k  =  fc(xn+1,xn+1)  (11) 


3..  Parameterizing  the  kernel  function 

There  are  many  possible  choices  of  prior  interpolating  kernel  functions.  From  a  modeling  point  of  view, 
the  objective  is  to  specify  a  prior  kernel  function  that  contains  our  assumptions  about  the  structure  of  the 
process  being  modeled.  Formally,  it  is  required  to  specify  a  function  that  will  generate  a  positive  definite 
kernel  matrix  for  any  set  of  inputs.  In  this  paper,  a  multi  layer  perception  (MLP)18  based  kernel  function 
is  used,  which  has  the  following  form 


h/  (xij  Xy  >  0)  —  O^Sin 


XT0>,+Qn 


VWWi  +  ^  +  i)(?$0%Xn  +  ^  +  1) 


+  *£ 


(12) 


In  Eq.  (12)  the  superscript  n  represents  the  nth  damage  instances.  It  is  noted  that  unlike  the  fixed  hyperpa¬ 
rameters,  the  hyperparameters  in  Eq.  (12)  are  adaptivly  found  as  new  data  set  D  =  {x^a*}^ . 

becomes  available. 


4>  Hyperpammeters  determination 

So  far  we  have  only  considered  the  properties  of  the  prediction  model  for  fixed  values  of  the  hyperparameters. 
This  section  discusses  how  to  obtain  the  hyperparameters  0  for  a  fixed  training  data  set  D  ■=  {xi,  a;}7=1  n. 
Ideally  integration  over  all  possible  hyperpajranieters  should  be  done  in  order  to  obtain  the  best  possible 
predictions  of  the  function  value  an+1  at  damage  level  n-f  L  Therefore,  it  can  be  written  as, 

f  (an+l\D,Xn+l,K(-))  =  J  f  (an+1|Ax„+i,  K (•),©)  /  (©|D, /<(<))<*©  (13) 

The  above  integral  is  as  complex  as  the  integral  given  in  Eq,  (6)  and  also  difficult  to  evaluate  for  a  complex 
problem  with  several  hyperparameters  and  a  multiple  input  space.  Out  of  the  two  possible  approaches  e.g., 
the  Maximum  evidence16  and  the  Moiite  Carlo17  approach  only  the  use  of  the  maximum  evidence  approach 
will  be  discussed  to  evaluate  the  integral.  Using  maximum  evidence  approach,  Eq.  (13)  can  be  written  in 
its  approximate  form  as 


/  K+i|£,xn-H,  K{-))  =  f  {ait+xlD-Xji+uK^.GMAp)  (14) 


The  approximation  in  Eq.  (14)  is  based  on  the  assumption  that  the  posterior  distribution  over  0,  he 
/  (0| £>,  ih(-))?  has  a  sharp  peak  around  Qmap  ♦  This  approximation  is  generally  reasonable13  and  predictions 
are  often  found  very  close  to  those  obtained  using  the  true  predictive  distribution.  Now  to  find  the  peak 
location  of  /  (S\D, JT(-))  the  posterior  of  it  needs  to  be  optimized  and  the  posterior  distribution  can  be 


written  as 


(15) 


Ill  Eq.(15),  the  denominator  (i.e.,  the  evidence)  is  independent  of  0  and  can  be  ignored  in  the  optimization 
process.  Oil  the  other  hand,  the  other  two  terms,  the  likelihood  /( an\  {x^}i ,  A"(*)?0),  and  the  prior  /  (0), 
need  to  be  considered  in  the  optimization  of  /  (0|D,  With  the  assumption  that  all  damage 

indices  follow  a  Gaussian  distribution  and  using  Eq. (8),  the  logarithm  of  the  objective  function  can  be  written 
as 


Log  {f(W>  £(')))=* 


| Log(dctKn ) 

\“ZKnlftn  ~  ^Log[2ix)  +  Logf  (©) 


(16) 
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The  log-likelihood  function  L  in  Eq.(16)  is  generally  multi-modal  and  can  be  optimized  using  any  multi¬ 
variate  optimization  algorithm.  In  the  present  work,  the  conjugate  gradient  method  is  used  to  optimize  the 
log-likelihood  function  and  to  obtain  the  optimized  hyp erparameters.  Note  that  it  is  common  practice13  to 
ignore  the  log  prior  term  in  Eq.(16)  due  to  the  absence  of  knowledge  on  0.  The  resulting  solution  may  not 
be  always  a  realistic  solution,  however  it  ean  be  assumed  that  Log f  (0)  is  implicitly  modeled  through  the 
optimization  of  the  log-likelihood  L . 

5.  Input- output  data  set  for  single  step  ahead  prediction 

For  single  step  ahead  prediction,  the  Gaussian  process  prediction  model  given  by  Eq.(8),  only  predicts  single 
step  ahead  damage  index.  For  prediction  of  the  n.+  1th  damage  index,  the  training  data  set  D  and  test  input 
vector  xn+i  can  be  stated  as, 


Training  data  matrix 


Target  vector 


ao 

a1 

0,2 

a>d 

ai 

«2 

ad 

ad+i 

an-d+l 

ari-d+ 2 

ton 

— 


T est  input  data  vector 


dfi  1  *  ■  •  &n 


(18) 


where  in  Eq.  (17  and  18)  the  subscript  n  symbolizes  the  nth  damage  instance  or  damage  level,  up  to  which 
the  last  on-line  data  was  available,  and  d  symbolizes  dimension  of  the  input  space. 


6.  Input-output  data  set  for  multi  step  ahead  prediction 

For  multi  step  ahead  prediction,  the  GP  model  given  by  Eq.(8)  recursively  predicts  the  future  state  multi 
step  ahead  of  the  last  on-line  data  available.  However,  unlike  the  single  step  ahead  prediction  case,  the  multi 
step  ahead  training  data  set  D  and  test  input  vector  xn±a  are  adaptively  updated  with  off-line  predicted 
damage  indices  rather  than  on-line  estimated  damage  indices.  For  prediction  of  the  n  4*  hth  damage  index 
an+n  =?  the  training  data  set  D  and  test  input  vector  can  be  written  as 


D 


[  I  ai  1 . 

L  J  i~a I4*n 


Training 

data  matrix 

Target  vector 

' 

a  o 

di 

a*2 

&d- 1 

Q>d 

ai 

d2 

a3 

*  &d 

&d- fl 

an-d 

tin-d+L 

&n—d+2 

* 

an 

ttn-d- j-1 

Ciji^d+2 

d‘Tt—d+3 

<+i 

a-n-d+ 2 

dn~d+3 

&n— d-f-4 

V 

■  <+l 

an+2 

P 

^TL—d—l+n 

ClP  . ,  ~ 
n—d±n 

ap 

d-fn+3 

np 

2+n 

an-l+n 

(19) 


^71 -Hi  — 


Test  input  data  vector 

.  . 

nP  np  np  np 

i  ^71— 2  *  *  *  •  n— 1+n 


(20) 
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where  in  Eq.  (19  and  20)  the  subscript  n  symbolizes  the  damage  instance  up  to  which  the  last  on-line  data  is 
available,  and  the  subscript  ft  symbolizes  the  damage  instance  number  following  the  availability  of  the  last 
on-line  data  available,  and  the  superscript  p  symbolizes  predicted  damage  index  from  the  off-line  module,  as 
opposed  to  being  estimated  from  the  on-line  model. 

7.  Residual  useful  life  estimation  (RULE) 

The  residual  useful  life  estimation  can  be  defined  as  the  difference  between  the  number  of  fatigue  cycles  at 
which  the  predicted  damage  index  becomes  critical,  Le  reaches  its  critical  value  (a*),  and  the  number  of 
fatigue  cycles  at  which  the  last  on-line  data  is  available*  The  RULE  can  be  defined  as: 

RULE  =  (n  +  n* )  AN  -  (n)  AN  -  N0  -  hAN  -  N0  (21) 

where  AN  is  the  number  of  fatigue  cycle  increments  per  each  increment  of  damage  instance  and  n*  corre¬ 
sponds  to  the  number  of  damage  instances  for  the  damage  index  to  become  critical  after  the  last  available 
on-line  data. 

Ill-  Numerical  Results 

A.  Fatigue  experiment  and  data  collection 

To  numerically  validate  the  integrated  prognosis  algorithm,  a  fatigue  test  was  performed  on  a  Al-6061. 
cruciform  specimen  under  biaxial  loading.  The  loaded  cruciform  specimen  in  a  MTS  biaxial  fatigue  test 
frame  can  be  seen  in  Fig.  3.  The  specimen  was  subjected  to  a  constant  amplitude  fatigue  loading  with 
maximum  amplitude  {crmQX)  4  kips  and  load  ratio  R— 0.1,  and  the  biaxial  machine  actuator  was  operated 
with  a  frequency  of  10  Hz.  It  is  noted  that,  the  maximum  stress  amplitude  omilx  —  4kips  is  equal  to  two 
third  of  the  yield  stress  ay.  Based  on  nonlinear  finite  element  analysis  of  cruciform  specimen,  the  yield 
stress  is  approximated  as  oy  —  Skips.  Also  note  that  both  the  x-axis  actuator  and  y-axis  actuator  of  the 
biaxial  frame  was  subjected  to  in-phase  fatigue  loading.  For  on-line  state  estimation,  passive  strain  gauge 
sensors  were  used.  One  strain  gauge  is  mounted  on  the  horizontal  flange  (Fig.  4a),  and  one  on  the  vertical 
flange  (Fig.  4a)  and  two  strain  gauges  are  mounted  on  the  web  area  (Fig.  4b)  of  the  cruciform:  specimen. 
In  addition,  to  the  strain  gauge  sensors,  piezoelectric  (PZT-5H)  sensors  were  also  instrumented  to  collect 
the  active  sensing  data.  However,  the  active  sensing  data  will  be  used  in  a  future  study.  The  current  paper 
is  based  on  strain  gauge  measurements.  A  hole  in  the  center  of  the  specimen  was  made  to  create  crack 
initiation  in  the  web  area  of  cruciform  specimen.  To  accelerate  damage  growth  an  EDM  notch  of  1  mm 
length  was  made  at  left  quadrant  boundary  of  the  central  hole  (45°  to  the  vertical  axis).  A  48  channel  NI 
PXI  system  was  used  to  collect  the  strain  gauge  signals,  the  piezoelectric  sensor  signals  and  the  measurements 
from  the  biaxial  machine  load  cells.  In  addition,  a  high  resolution  SONY  camera  was  also  used  to  visually 
monitor  the  crack  growth.  The  data  acquisition  system  and  the  computer  capturing  the  visual  image  were 
synchronized  with  the  biaxial  machine  controller  to  collect  the  time  synchronized  data/  image  at  a  specified 
interval  of  AN  =  1500c?/cles.  The  data  and  image  collection  started  at  approximately  Ilk  cycles.  The 
image  and  sensor  data  was  collected  at  47  different  time  instance.  For  the  first  44  instances,  the  signals  and 
images  were  collected  while  the  biaxial  machine  was  running  and  during  the  last  three  instances  the  data 
was  collected  when  the  machine  was  not  running.  This  lead  to  a  total  of  44  different  damage  cases  with  the 
last  damage  state  occurring  at  75.5  kcycles.  The  developed  real-time  MATLAB  based  prognosis  algorithm 
was  also  synchronized  with  the  data  acquisition  system  to  estimate  the  current  damage  state,  and  to  predict 
the  future  damage  state  and  remaining  useful  life. 

B.  Correlation  features  and.  damage  index 

To  evaluate  the  damage  index  as  mentioned  earlier,  strain  signal  are  mapped  as  input  and  output.  For 
example,  the  signal  (e^)  from  the  strain  gauge  mounted  on  the  horizontal  flange  (or  X-arm)  of  the  cruciform 
specimen  is  considered  as  the  input  signal  u,  whereas  the  signal  (ej^ )  from  the  web  mounted  stain  gauge  is 
considered  as  output  v.  It  is  noted  that  both  the  horizontal  axis  strain  (e^)  and  the  vertical  axis  strain  (ej^) 
are  measured  by  two  different  strain  gauges  placed  perpendicular  to  each  other  .  Comparison  of  input  strain 
(e£)  and  output  strain  (e]f)  is  shown  in  Fig.  5.  It  presents  a  comparison  for  four  different  damage  cases, 
damage  case  7  (at  20  kcycle),  damage  case  20  (at  39.5  kcycle),  damage  case  42  (at  72.5  kcycle)  and  damage 
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Figure  3.  Al-6061  cruciform  specimen  loaded  in  a  MTS  biaxial  fatigue  test  frame 


case  44  (at  75.5  kcycle).  From  the  figure  it  can  be  seen  that  though  there  is  a  clear  trend  between  input 
and  output  strain  during  the  later  stages,  the  trend  is  hardly  discernabie  in  the  earlier  stages  of  damage 
growth.  Rather  than  directly  using  the  time  series  data  for  different  damage  case  comparisons,  using  Eq. 
(4),  the  cross-correlation  coefficient  between  input  and  output  is  found  for  different  damage  cases.  The 
comparison  of  cross-correlation  coefficients  for  damage  level  1  (reference  case  at  11  k  cycles)  with  cross¬ 
correlation  coefficient  at  different  damage  levels  are  shown  in  Fig.  6.  Figure  6a,  6b,  6c,  and  6d,  respectively 
show  the  comparison  of  cross-correlation  coefficients  of  damage  case  1  with  damage  case  7,  damage  case  20, 
damage  case-42  and  damage  case-44.  It  is  noted,  the  results  shown  in  Fig.6,  the  x-axis  flange  strain  (ef ) 
and  x-axis  web  strain  e ™ )  are  respectively  taken  as  input  u  and  output  v.  Also,  it  can  be  seen  there  is  a 
better  trend  of  damage  condition,  compared  to  the  direct  time  scries  data  from  the  strain  gauges.  However, 
to  compare  the  different  damage  levels  using  a  scalar  quantity  and  to  improve  the  discernability  between 
different  damage  cases,  the  damage  index  proposed  in  Eq.  (5),  is  evaluated  for  the  different  damage  states. 
Figure  7  shows  increase  in  the  damage  index  with  fatigue  cycle  for  two  different  output  measurements, 
and  e™  y  for  fixed  input  from  x-axis  flange  measurement  ef .  The  figure  shows  a  clear  trend  of  damage 
growth  with  e ™  as  output  strain.  A  good  correlation  between  estimated  damage  index  and  normalized  visual 
measurement  is  also  observed.  It  must  be  noted  that,  the  visual  measurement  is  available  up  to  damage  level 
29  (up  to  53  kcycle).  After  the  29th  damage  level,  it  was  found  that  the  camera  went  out  of  focus.  Figure 
7  also  shows  that  with  respect  to  as  the  output  strain,  except  for  inal  failure  regime,  there  is  no  clear 
trend  in  damage  growth.  This  is  because  the  input  signal  c£  is  poorly  correlated  with  the  y-axis  web  strain 
(cjn  measurements.  The  higher  damage  indices  during  the  final  failure  regime  is  possibly  due  to  presence 
of  shear  strain  components.  In  addition  to  the  strain  signal  as  input,  a  damage  index  is  also  obtained  using 
biaxial  frame  load  cell  measurements  (x-axis  load  cell).  The  new  damage  index  is  shown  in  Fig.  8,  and  a 
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Figure  4.  Undamaged  and  damaged  condition  of  cruciform  specimen:  Plot-a  shows  the  undamaged  cruciform 
specimen.  This  rear  view  of  the  specimen  also  shows  the  location  of  two  strain  gauges  mounted  in  the  web 
area.  Plot-b  shows  the  final  damage  condition  (at  75.5  kcycles)  of  the  cruciform  specimen.  This  front  view 
of  the  specimen  also  shows  the  location  of  two  strain  gauges:  one  mounted  on  horizontal  arm  and  the  other 
mounted  on  the  vertical  arm  of  the  specimen. 


similar  trend  in  damage  index  growth,  as  in  the  case  of  x-axis  flange  strain  measurement  as  input,  can  be 
observed.  However,  it  is  noted  that  in  a  real  life  system,  it  is  hardly  possible  to  directly  measure  the  loads 
applied  to  the  structure.  On  the  other  hand,  it  is  feasible  to  mount  strain  gauges  or  small  sensors  at  required 
locations  without  affecting  the  structural  integrity  of  the  host  structure.  Therefore  the  results  presented  in 
the  subsequent  sections  are  based  only  on  the  strain  gauge  based  data. 

C.  Single  step  ahead  state  prediction 

Figure  9  shows  the  comparison  between  single  step  ahead  prediction  and  actual  damage  index  with  on¬ 
line  data  available  up  to  the  previous  damage  level.  As  seen  in  the  figure,  the  prognosis  algorithm  starts 
predicting  from  damage  level  7.  It  is  noted  that  the  dimension  d  of  the  Gaussian  process  input  is  chosen 
as  6.  Therefore  the  prognosis  algorithm  requires  at  least  six  damage  states  to  obtain  the  1x6  test  input 
vector  (see  Eq.  18).  Also,  with  unavailability  of  any  training  data  set  D  (Eq.  17),  to  predict  the  7th  damage 
state,  the  initial  hyperparameters  (Eq.  16)  are  chosen  as:  =  0™  =  1  and  0%  =  0.1.  Because  of  this, 

there  is  a  large  mismatch  between  the  7th  level  damage  prediction  and  the  actual  damage  index.  However, 
for  prediction  of  damage  level  eight  and  beyond,  the  traing  input  data  matrix  (Eq.  17)  and  target  vectors 
(Eq.  18)  arc  recursively  updated.  For  each  recursive  updating,  a  new  set  of  hyperparameters  arc  obtained 
using  the  conjugate  gradient  method.  Once  the  hyperparameters  are  estimated,  the  one-step  ahead  damage 
index  is  predicted  for  the  following  damage  level.  Figure  9  shows  a  clear  correlation  between  one-step  ahead 
predicted  damage  level  and  the  actual  damage  index  before  the  damage  index  reaches  its  critical  value  of  0.7. 
It  must  be  noted  that  the  actual  damage  indices  arc  the  on-line  damage  states  (or  damage  index);  these  are 
directly  estimated  from  the  sensor  signals.  The  threshold  value  of  0.7  is  70%  of  the  final  damage  index  value 
of  1.  From  Eq.  (5)  the  damage  index  reaches  its  final  value  of  1  when  there  is  no  cross-correlation  between 
the  input  u  and  v.  This  is  because  the  specimen  undergone  complete  failure.  It  is  noted  that  choosing  the 
critical  damage  index  value  of  0.7  is  based  on  the  results  from  previously  performed  similar  experiments. 

D.  Multi  step  ahead  prediction 

Unlike  the  single  step  ahead  prediction,  the  multi  step  ahead  prediction  recursively  predicts  the  damage 
state  ahead  of  the  last  online  data  availability  stage.  Figure  10  shows  the  multi  step  ahead  state  prediction. 
Similar  to  single  step  ahead  prediction  process,  the  prognosis  algorithm  initiates  after  the  6th  damage  level 
(at  18.5  kcycles).  From  the  7th  damage  level  (from  20kcyclcs),  damage  indices  are  predicted  and  then  fed 
back  to  the  prognosis  model  to  update  the  training  data  matrix  (Eq.  19)  and  the  test  input  vector  (Eq.  20). 
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Figure  5.  Input  output  strain  comparisons  at  different  damage  levels. 


The  feedback  process  is  continued  recursively  as  long  as  the  predicted  damage  index  does  not  reach  its 
critical  value  of  70%  of  its  final  failure  value.  It  is  to  be  noted  that  unlike  single  step  ahead  prediction,  the 
training  data  matrix  and  the  test  input  vector  are  updated  with  off-line  model  predicted  states,  rather  than 
being  updated  with  on-line  model  estimated  states.  It  can  be  seen  from  Fig.  10  that  up  to  damage  level 
23  (at  44  kcycles)  the  multi  step  ahead  prediction  fails  to  reach  the  critical  damage  index.  This  is  because 
the  predictive  model,  which  is  unable  to  learn  the  damage  growth  dynamics.  To  reduce  the  computational 
expenses,  the  prognosis  algorithm  was  stopped  at  certain  times.  The  criteria  for  stopping  the  algorithm  was 
if  the  rate  of  damage  index  growth  is  not  greater  than  1  x  10 ~7 /cycles  for  six  consecutive  damage  levels, 
the  off-line  predictive  model  is  terminated.  This  is  because,  physically,  if  the  damage  growth  is  so  slow  the 
damage  index  value,  will  never  reach  the  critical  value  even  if  the  algorithm  were  to  run  indefinitely.  From 
Figure  10  it  is  also  seen  that,  the  first  multi  step  ahead  prediction  curve,  that  reaches  the  critical  value  starts 
from  damage  level  24  (from  45.5  kcycles).  Beyond  this  damage  level,  the  multiple  step  ahead  prediction 
increasingly  converges  with  the  actual  damage  index.  Also,  it  can  be  assumed  that  the  prediction  horizon 
is  between  damage  level  24  (45.5  kcycles)  and  damage  level  42  (72.5  kcycles),  at  which  the  on-line  model 
estimation  states  reaches  its  critical  value. 

E.  Residual  useful  life  estimation  (RULE)  and  mean  square  error  evaluation 

Using  Eq.  (21),  the  residual  useful  life  at  a  given  damage  level  (up  to  which  the  last  online  data  was  available) 
is  estimated.  Figure  11  shows  the  comparison  of  predicted  RULE  and  actual  RULE.  From  the  figure  it  can 
be  seen  that,  there  is  a  good  correlation  between  predicted  and  actual  RULE  in  the  true  positive  regime  i.e., 
between  damage  level  24  (at  45.5  kcycles)  and  damage  level  42  (at  72.5  kcycles).  Also  it  can  be  observed 
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Figure  6.  Cross-correlation  coefficient  comparison  at  different  damage  levels  with  the  reference  level  at  11 
kcycles 


Figure  7.  Variation  of  damage  index  with  fatigue  cycle.  Flange  (x-axis)  strain  measurements  are  used  as  input 
and  web  (x-axis)  strain  measurements  as  output 
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Figure  8.  Damage  index  with  x-axis  biaxial  frame  load-cell  measurements 
measurements  as  output 


as  input  and  web 


(x-axis)  strain 


that  as  more  and  more  online  data  becomes  available,  better  correlation  between  predicted  RULE  and  actual 
RULE  is  obtained.  Figure  12  shows  the  mean  square  error  between  predicted  RULE  and  estimated  RULE. 
It  can  be  seen  that  during  the  true  positive  regime,  the  mean  square  error  is  substantially  reduced  compared 
to  the  mean  square  error  during  the  false  positive  regime. 
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Figure  9.  Off-line  model  prediction  of  one-step  ahead  damage  state. 


IV .  Conclusion 


An  adaptive  on-line  off-line  life  integrated  predictive  model  has  been  developed  for  applications  to  Struc¬ 
tural  Health  Monitoring  and  damage  prognosis.  The  developed  prognosis  model  combines  an  on-line  state 
estimation  model  with  an  off-line  predictive  model  to  adaptively  estimate  the  residual  useful  life  of  an  Al- 

cruciform  specimen  under  biaxial  loading.  Some  important  features  of  this  algorithm  and  observation 
irom  this  study  arc  summarized  below. 


1.  The  on-line  model  is  based  on  a  non-parametric  system  identification  approach,  which  estimates  the 
cun-ent  damage  states.  The  numerical  results  showed  good  correlation  between  on-line  estimated  state 
and  the  normalized  visual  measurements. 

2.  Once  the  current  damage  state  was  available  from  the  on-line  model,  the  information  was  fed  to  an 
o  -line  predictive  model  to  obtain  the  future  states  and  remaining  useful  life  estimation  (RULE). 

3.  The  off-line  predictive  model  is  a  high-dimensional  kernel  function  based  recursive  Gaussian  process 


4.  The  future  states  are  recursively  predicted  by  feeding  back  the  previous  predicted  states  to  the  off-line 
model.  Also,  the  model  parameters  (Gaussian  process  hyperparameters)  are  updated  with  repetitive 
conjugate  gradient  based  optimization. 


5.  Good  correlation  was  observed  between  actual  damage  states  and  predicted  future  damage  states 
between  the  24  damage  level  (at  45.5  kcycles)  and  the  41st  (at  71  kcycles)  damage  level  at  which  the 
damage  index  reaches  its  critical  value  of  0.7. 


6.  The  duration  between  the  24th  damage  level  (at  45.5  kcycles)  and  the  42"d  damage  level  (at  72  5 
kcycles)  was  considered  as  true  positive  regime  or  prediction  horizon.  This  regime  also  showed  that  a 
good  correlation  between  predicted  RULE  and  actual  RULE. 
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Figure  10.  Off-line  model  prediction  of  multi  step  ahead  damage  state. 
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Figure  11.  Comparison  of  predicted  RULE  and  actual  RULE. 
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Figure  12.  Mean  square  error  between  predicted  RULE  and  actual  RULE. 
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